Bitwise Bell inequality violations for an entangled state involving 2^ ions 
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Following on from previous work [J.-A. Larsson, Phys. Rev. A 67, 022108 (2003)], Bell inequalities 
based on correlations between binary digits are considered for a particular entangled state involving 
2N trapped ions. These inequalities involve applying displacement operations to half of the ions 
and then measuring correlations between pairs of corresponding bits in the binary representations 
of the number of centre-of-mass phonons of N particular ions. It is shown that the state violates 
the inequalities and thus displays nonclassical correlations. It is also demonstrated that it violates 
a Bell inequality when the displacements are replaced by squeezing operations. 
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I. INTRODUCTION 

Entangled quantum states typically, if not always, ex- 
hibit nonclassical correlations. These correlations are 
crucial elements in most quantum information process- 
ing tasks including quantum computation 0,0, quantum 
teleportation Q, superdense coding [4| and some forms 
of quantum cryptography [ij (Section 12.6). Given this 
significance, it is important to consider how to best ob- 
serve such correlations and thus better understand the 
quantum resources present in certain situations. One 
way of observing nonclassical correlations is via the viola- 
tion of Bell inequalities [MEl- For example, violations of 
the Clauser-Holt-Shimony-Hornc (CHSH) inequality Q 
can reveal the presence of such correlations in many two- 
qubit entangled states. Similarly, a violation of the GHZ 
inequality || || highlights the existence of nonclassical 
correlations in the state IV'ghz) = 1/V2 (|000) + | 111)). 
Finally, the violation of a Bell inequality involving higher- 
dimensional spin [13 by the spin-s singlet state, where 
s = 3/2,5/2,7/2... highlights nonclassical correlations 
present in this state. Many other Bell violations are also 
known, however, they are too numerous to mention. 

The examples in the previous paragraph involved Bell 
inequalities well-suited to observing nonclassical correla- 
tions in particular entangled states. However, not all 
Bell inequalities are useful for observing such correla- 
tions in every entangled state. For instance, apply- 
ing the CHSH inequality to any two qubits in IV'ghz) 
produces no Bell violation and hence, when used in 
this manner, this inequality does not highlight It/'gh zVs 
nonclassical correlations. Similarly, the W state [13 
|V>w) = l/v / 3(|001) + |010) + |01)) satisifies the GHZ 
inequality and hence this inequality is ill-suited for ob- 
serving its nonclassical correlations. Another noteworthy 
point about Bell violations and entangled states is that 



certain mixed entangled states, namely bound entangled 
states, may not violate any Bell inequality as has been 
conjectured by Peres Consistent with this, it has 

been shown that multipartite bound entangled states for 
which all partial transposes are positive satisfy one par- 
ticular Bell inequality |12| . 

Given that individual Bell inequalities can be either 
good or bad tools for observing nonclassical correlations 
in specific entangled states, it seems interesting to con- 
sider the following question: "Which particular Bell in- 
equalities are best suited to observing the nonclassical 
correlations of a certain entangled state?" Whilst not 
addressing this general question in the current paper, we 
do show that certain Bell inequalities involving correla- 
tions between binary digits in the binary representations 
of particular observables can be used to detect interesting 
nonclassical correlations in a particular entangled state 
involving two sets of N ions. In doing so, we follow on 
from [13| which showed that the steady-state intracavity 
state of the nondegenerate parametric amplifier (NOPA) 
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where \n) is a photon number state, the subscripts 1 
and 2 denote the signal and idler modes and r is a 
real squeezing parameter, violated certain Bell inequal- 
ities. More specifically, using an abstract mathematical 
scheme 0] showed that if we consider the numbers of 
photons in the signal and idler modes in binary (eg. 
I ra = 3)i 1— > 1 71 = ...00011)i) then each pair of corre- 
sponding bits in the two binary representations simul- 
taneously violates a CHSH Bell inequality. That is, it 
showed that the least significant bits for the signal and 
idler modes, together, violated such an inequality as did 
the second least significant bits, the third least significant 
bits and so forth. The paper [lj| also briefly suggested 
how we might observe these violations but, on this point, 
remarked that a better (that is, presumably, a more ex- 
perimentally achievable) measurement scheme than that 
suggested was desirable 13] (p. 022108). 
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The results in [13| can be seen as extending those in 
El 13 El 13 El which all showed that |NOPA) violated 
Bell inequalities based on measuring photon-number par- 
ity (oddness or even-ness) . In terms of binary represen- 
tations, these other papers violated Bell inequalities in- 
volving the values of the least significant bits (and not 
any other bits as did [13|) m the binary representations 
of the numbers of photons in the signal and idler modes. 

This paper extends and complements work in EH by 
explicitly showing the existence of Bell violations closely 
related to those in |13| within a tangible and arguably ex- 
perimentally feasible context that is different to the con- 
text suggested in EH • In particular, motivated by EH ' s 
comment (p. 022108) that the formulation of a practi- 
cal measurement scheme for the Bell inequalities in EH 
is desirable, we (arguably) propose such a scheme. This 
paper also extends work in [TJ| by illustrating different 
ways to violate the sorts of Bell inequalities in to the 
ways shown in |13( . Instead of violating Bell inequalities 
by measuring a range of pseudo-spin observables for the 
state |NOPA) as did EH> we a PPby a range of displace- 
ments and squeezing operations to a state generated from 
|NOPA) and then always measure the same pseudo-spin 
observables. 

The current paper proceeds as follows: in Section [H] 
the state we consider is described, along with the phys- 
ical system underlying it which centres around two sets 
of N trapped ions. In Section IIIII the Bell inequalities 
we consider are presented by outlining the measurements 
and operations they involve. The measurements consist 
of measuring bits in the binary representations of A/"i and 
A/2, where A/j (j = 1,2) is the number of centre-of-mass 
phonons for one of the sets of N ions in the x direc- 
tion, whilst the operations (which are applied prior to the 
measurements) are displacements applied to the centre- 
of-mass vibrational states of one of the sets of ions. In 
Section ITV1 it is shown that the entangled state violates 
the Bell inequalities. Next, Section El presents a Bell in- 
equality involving local squeezing operations which the 
entangled state also violates. Finally, the paper con- 
cludes with a discussion of its results in Section IVT1 



II. THE ENTANGLED STATE 

The system associated with the entangled state con- 
sidered in this paper comprises of a nondegenerate para- 
metric amplifier (NOPA) [3 El EH and two linear ion 
traps which each lie within an optical cavity and contain 
N identical ions. A schematic diagram of this system 
is shown in Fig. ^ The NOPA operates below thresh- 
old and its two external output fields first pass through 
Faraday isolators. Each of them then feeds into a dif- 
ferent linearly-damped optical cavity via a lossy mirror. 
The cavities are aligned such that their axes coincide with 
the x axis and are closed at one end by perfectly reflect- 
ing mirrors. In addition, each cavity supports a cavity 
mode of frequency lo c described by the annihilation op- 



erator Oj, where j enumerates the cavities. Within both 
cavities lie N identical two-level ions of mass M, charge 
Z and internal transition frequency oj a . These ions are 
trapped in a linear configuration parallel to the x-axis by 
a harmonic potential (a linear ion trap and hence 
are tightly confined in the y and z directions. Further- 
more, the vibrational motion of the m th ion in the j 
trap in the x direction is described by the annihilation 
operator for which ^] = 1. The traps are 

aligned such that the j th trap is centred on a node of 
the cavity field described by aj. Finally, external lasers 
of frequency ujl whose beams are perpendicular to the 
x-axis are incident on the first ions of both traps. 
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FIG. 1: Schematic diagram for the system associated with 
the 2iV-way entangled state. The system consists of, firstly, 
a subthreshold optical nondegenerate parametric amplifier 
(NOPA) whose output modes pass through Faraday isolators 
(represented by an F enclosed in a circle) and then feed into 
linearly damped optical cavities. These cavities are aligned 
along the a;- axis and each has one ideal mirror and one lossy 
one (with damping constant n). Inside each cavity is a har- 
monic ion trap that confines N identical two-level ions (each 
represented by a black circle) in a linear chain parallel to the 
x-axis. External lasers of frequency ujl are incident on the 
first ions in both traps from a direction perpendicular to the 
x-axis. 

The Hamiltonian for the j th optical cavity, the ions 
within it and its reservoir is 

tt Tjion 1 [f 1 frion— ion , tt 

-Hj total — tljo + -"j'O + -Hjj + ti.ji 

+n{ aj R] + a]Rj) + H } res , (2) 

where H^ n is the free Hamiltonian for the vibrational 
states of the ions and Hjo is the free Hamiltonian for 
the cavity field and the ions' internal states. The 
term £fj° n-lon describes the electromagnetic coupling be- 
tween ions whilst Hjj describes a Raman process in- 
volving the cavity field, the external laser and the first 
ion in the j th trap. Finally, Hj res is the Hamilto- 
nian for the external reservoir coupled to the j th cav- 
ity for which Rj is a reservoir annihilation operator 
and k is a damping constant. More precisely, = 

ELi (^r ) t& j™ } + 5) > where v x (v x = uj c - uj l ) is 
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the frequency of both harmonic traps along the x-axis. 
The Hamiltonian Hj is, in a frame rotating at frequency 

N 

H j0 = hSat aj + hAj2 ojf o-}- 5 > (3) 

m— 1 

where 5 = uj c — u>l, A = uj a — ujl, and and are 
raising and lower operators for the internal states of the 
m th ion in the j th trap. The term -ffj-° n-ion is 

7 2 

rrion— ion y^N f ax 

jl m,n— 1: m^n o I /i \ i ' V / 

where eo is the permittivity of free space and Xji, for 
I = 1 ... N, is the position of the I th ion in the j th trap. 
The interaction term Hji is 

Hjj = H[£ L (y,z,t)a^+£* L (y,z,t)af2] 

+hg sin{kx J i)(a]cjf} + ajofj), (5) 

where is the complex amplitude for both external 
lasers, k — lo c /c and go ( go ^ 5ft) is the coupling constant 
for the ion-field interaction. 

The following feasible assumptions are made about the 
system |2(| in order to simplify calculations and to focus 
on its most important aspects: 

1. All ions are so cold that they only move from their 
mean position x^ by a small amount and so we can 
approximate Xji(t) by x^ + Qji(t), where qji(t) is a 
small displacement. 

2. The cavity field and external laser frequencies are 
appreciably detuned from u a and all two-level ions 
are initially in their ground states. Thus, the ex- 
cited internal states are sparsely populated and 
spontaneous emission effects are negligible and can 
be ignored. 

3. The wavelength of the cavity mode is much greater 
than the distances that the first ions in both traps 
stray from the centres of their traps and thus 
sin(fcxji) ~ kxji << 1. This allows us to to ar- 
range things so that the y and z dependences of 
the external laser fields are negligible and thus, as- 
suming £l is time independent, that £L(y,z,t) ~ 
£ e -40L ; where £ is a real time-independent ampli- 
tude. 

4. The damping parameter k is such that v x ^> k 
g k£Vh/ (y/2Mv x NA) . 

5. For each trap, the frequencies of different normal 
or collective modes [2~H in the x-direction are well- 
separated. Thus, the cavity modes only couple to 
the centre-of-mass modes in this direction. 



Given assumptions 1, 3 and 5, calculations in |25| show 
that we can write Hj total m terms of normal-mode cre- 
ation and annihilation operators as 

Hi tota, = fcf> m (b]MbM + 1)+ H J0 
m=l ^ ' 

+h[£ L (y,z,t)<j^+£* L (y,z,t)^ 1 }] 

+K(ajRt + a]Rj) + H } ros , (6) 

where is the annihilation operator for the m th nor- 
mal mode for the j th trap in the x direction. For exam- 
ple, Bj x is a centre-of-mass mode annihilation operator 

which is Bf] = 1/VN (b? x ] + bf x +... 6^) whilst Bf] 

is the annihilation operator for the breathing mode which 

is Bf x ] = 1/V2 (-bf x ] + bfj^j when N = 2. Observe that 

in Eq. © the cavity mode only couples to the centre-of- 
mass vibrational mode in the x direction. 

Though it would be very challenging, at best, to ex- 
perimentally realise the system outlined above, it is po- 
tentially feasible to do so. This is because, first, optical 
cavities and parametric oscillators have been widely re- 
alised in laboratories. Second, recent experiments |27j 
have trapped single ions in electromagnetic traps lying 
within optical cavities. 

Upon adiabatically eliminating the cavity mode a,j and 
also cr,-_|_ and (jj- in Eq. iJBJ, it can be shown that the 

Langevin equation for B^J is 

where x denotes the partial derivative of x with respect 
to time, r = hqnk 2 £ 2 / (2Mv x kA 2 ) and aj n is a quantum 
noise operator 30] . This equation shows that the only 
effect of having multiple ions, as opposed to a single ion, 
in the trap is to introduce a factor of 1/N in front of 
r. From [26j (especially Eqs (11) and (12)), it is known 
that when N = 1 the evolution described by Eq. (JTJ 
implements a process known as quantum state exchange 
|26j which involves the transferral of the quantum state 
of an electromagnetic field mode to that of one or more 
trapped atoms. In particular, it is known that, for N — 
1, Eq. J7J) implements this process via the transferral of 
information about the input field a- lTL to the ion's centre- 
of-mass vibrational state. It thus follows that Eq. Q 
also implements quantum state exchange when N > 1, 
albeit more slowly due to an effective decrease in F with 
N in Eq. 0. 

In |23, it was shown that we can transfer the intra- 
cavity steady-state for the subthreshold nondegenerate 
parametric amplifier |NOPA) into the vibrational states 
in the x direction for two single trapped atoms in differ- 
ent harmonic traps. Using the connection between the 
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quantum state exchange processes involving a single har- 
monically trapped ion and N harmonically trapped ions 
demonstrated above, it follows that for the system illus- 
trated in Fig. n we can transfer |NOPA) into the centre- 
of-mass modes in the x direction of the two sets of N 
trapped atoms thus producing, in the steady state, 
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where \AP)j denotes the centre- of-mass vibrational num- 
ber state for the x direction with eigenvalue Af for the 
ions in the j th trap. 

For the remainder of the paper, we consider IV'cm) and, 
in Subsection IIII Bl present certain Bell inequalities in- 
volving correlations between bits in the binary represen- 
tations of A/"i and Af 2 . Section Hvl then shows that IV'cm) 
violates these inequalities and thus that they highlight 
nonclassical correlations in IV'cm). 

We have assumed the trapped ions suffer no vibrational 
decoherence. This assumption is justified in the following 
sense: In realistic systems, the timescale over which ap- 
preciable vibrational decoherence occurs is greater than 
that over which quantum exchange would take place 
[2HE1- Because of this, we can, in principle, generate a 
state very similar to | Vcm) before appreciable vibrational 
decoherence has occurred and then consider this state. 
Hence, even taking vibrational decoherence into account, 
we can produce a state very close to |Vcm), thus allowing 
us to ignore this decoherence of the trapped ions in our 
analysis. 



III. SCHEME FOR BITWISE BELL 
INEQUALITIES 

A. Motivation 



Let us for the moment pretend that each bit in Eq. © 
represents a physical qubit, with corresponding bits in 
each statevector with a '1'('2') subscript representing the 
same qubit. That is, with all of the least significant bits 
in each statevector denoted by a '1' ('2') subscript repre- 
senting one qubit, the second least significant bits in each 
statevector denoted by a '1' ('2') subscript representing 
another and so forth. Upon adopting this fiction, we see 
that |Vcm) factorises as follows: 



IV' 
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CM/ 
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where the superscripts label pairs of qubits. Performing 
single-qubit rotations on all qubits in Eq. IjlUfl and then 
making measurements in the computational basis, we can 
concurrently violate the CHSH inequality for all qubit 
pairs denoted by the same superscript. That is, we can 
simultaneously violate this inequality for the qubit pair 
denoted by (0), the one denoted by (1) and so forth. 

To observe IV'cm) 's nonclassical correlations we would 
like to implement the scheme involving CHSH violations 
described in the previous paragraph as it produces the 
largest possible violation for each qubit pair. However, 
our physical system of interest does not have the distinct 
qubits used in the scheme and so, in practice, we cannot 
address all binary digits individually. In spite of this we 
can still implement a similar scheme using other local 
unitaries and other measurements to produce multiple, 
though smaller, CHSH violations, as shown in the next 
subsection. 



B. The scheme 



To ease the reader into the Bell inequalities we con- 
sider, we now outline a line of thinking by which some- 
one might come to consider the closely related bitwise 
Bell inequalities in |l3| . 

Each centre-of-mass vibrational number state consti- 
tuting IV'cm) can be expressed as an infinite- length bi- 
nary string that denotes the number of centre-of-mass 
phonons in the state 0|. For example, |A/i = 2)i can 
be written as \Afi = ...0010)i. Expressing all centre- 
of-mass vibrational number states in this manner, IV'cm) 
becomes 
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tanh 2 r|...010)i| 
tanh 3 r|...011)i|... 011)2 + ...) 



..000)2 
.001) 2 
..010)2 



(9) 



where it is implicit that the bit strings represent Af% or 
A/2 values. 



In this subsection we present three CHSH inequalities 
involving correlations between three pairs of bits in the 
binary representations of A/i and A/2 ■ These inequalities 
involve displacement operations which we apply to both 
sets of ions before making certain measurements involv- 
ing electronic states. 

Applying displacements to both sets of ions in IV'cm) 
yields 



IVd) = D 1 (a)D 2 ((3)S 12 (r)\0) 1 \0) 2 



(11) 



where Di(a) and D 2 ((3) are, respectively, displacement 
operators acting on the first and second sets of ions in 
IV'cm) (that is the sets in the first and second traps, re- 
spectively) with displacements a and (3. These opera- 
tors are given by Di(a) — exp(aB[ 1 J ^ — a*B[ x J) and 

D 2 (j3) = exp(/3B^ )f - 0*B^). The operator S 12 (r) 
is the two-mode squeezing operator which is, when r is 

real, S 12 (r) = exp (r(5« f - B^B^)) , where r 

is a squeezing parameter. Lastly |0) 1 10) 2 is the two-mode 
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vacuum state for the centre-of-mass modes of the first and 
second sets of ions in the x direction. 

After applying D\ (a) and D2 ((3) , the next step in our 
CHSH inequality violations is to measure the values of 
the N least significant bits of A/i and A/jj ■ This is done 
using the measurement scheme in |32j which we now de- 
scribe. This scheme measures the N least significant 
bits of the number of centre-of-mass phonons for a set 
of N identical two-level ions (with internal transition fre- 
quency loq) in a linear ion trap. It begins by first set- 
ting the state of each ion to be an equal superposition 
of its ground and excited internal states. Next, the mea- 
surement scheme involves applying a standing-wave laser 
pulse to each ion such that each ion's mean position co- 
incides with a node of its pulse. The laser frequency for 
all pulses ojl is far detuned from all resonant vibrational 
frequencies for the trapped ions and A', the detuning 
between ujl and ujq, is such that |A'| 3> v x , where v x 
is the trap frequency in the direction along which the 
ions are aligned. The m th laser pulse is applied to the 
corresponding ion for a time t m — 2 m TTNA'/(2 N rilQ 2 ), 
where r\ x is the Lamb-Dickie parameter common to all 
ions and f2 is the Rabi frequency for each ion. Finally, 
an inverse Fourier transformation is applied to the ions' 
internal states. The measurement scheme has the effect 
of transferring the value of the m th bit of the number of 
centre-of-mass phonons for the ions to the two-level in- 
ternal system of the m th ion. This bit is encoded using 
the following mapping : \g) m <— > and |e) m 1— * 1, where 
\g)rn and \e) m are the ground and excited internal states 
for the m th ion. Once encoded in internal states, the bit 
can be readily measured using the resonant fluorescent- 
shelving technique 33] . 

Applying the -Di(a) and Z^/?) to \iPd) and then us- 
ing the measurement scheme on both sets of ions in the 
resulting state produces 

hM = 

00. 00 

£ c *Mf h \M h ® |binary(z, 1)}^ . . . |binary(i, N))\ <N 

1, j=0 

®|binary(j, 1))^ . . . |binary(j, N))l >N , (12) 

where binary(x, y) is the value of the y th 
least significant binary digit of x and Ci_j = 
(i,j\D 1 (a)D 2 ((3)S 1 2(r)\0) 1 \0}2. Superscript e's and 
vib's denote, respectively, an internal (or electronic) 
state and a centre-of-mass vibrational one for the x 
direction. Finally, the subscript k and / in |)| l denote 
that the state is for the I th electron in the fc th set of ions. 

We now assume, for the moment, that all observable 
quantities in \iPf) behave classically and thus can be 
simulated using a local hidden variable theory (LHVT). 
Given this, it follows that the correlations between the 
y th least significant bits of A/i and A/2, where y — 
1, 2, 3 ... N, can be described by a LHVT theory. Hence, 
using reasoning in Q, these correlations satisfy the 



CHSH inequality: 

S y = \E[a y (a)(b y ((3)+b y (f3>)) (13) 
+a y {a'){b y {(3)-b y {(3'))} | < 2, 

where a y (z) and b y (z') are the values of the y th least sig- 
nificant bits of, respectively, the first and second sets of 
ions given either Di(zi) or 1)2(22), where Z\ — a, a' and 
Z2 = (3,(3'. The notation E[. . .] denotes an average or 
expectation value. As Inequality (|13fl involves thinking 
about A/"i and A/2 binary digit by binary digit, we call 
the inequality in this equation a bitwise Bell inequality. 
Inequality H13|) arises from the fact that LHVTs are com- 
mitted to the existence definite values for all a y and b y 
at all times that can only change in a local manner. 

One important feature about the Bell-inequality 
scheme outlined above is that it is potentially realis- 
tic. This is because, first, the application of Di(z\) and 
1)2(22) to IV'Cm) is feasible as existing experiments have 
applied such operations to the vibrational state of a single 
trapped ion (see, for example, H^])- Second, the scheme 
is potentially realistic as the interaction between internal 
and centre-of-mass vibrational states in the measurement 
scheme it requires seems to be experimentally feasible. 
This is the case as it only requires far-detuned stand- 
ing wave laser pulses that interact with a particular ion 
for set times. Finally, it is conceivably feasible as the 
resonant fluorescent-shelving technique it uses to make 
measurements on internal states has been experimentally 
implemented with high efficiency (see, for example, |33jp. 



IV. RESULTS 

In this section, we show that IV'Cm) violates the three 
bitwise Bell inequalities represented by Eq. (|13f) when 
y = 1,2 or 3. Thoughout, we assume that N > 3 and 
hence that the measurement scheme can measure up to, 
at least, the third least significant bits in the binary rep- 
resentations of A/i and A/2 ■ 



A. Least significant bits 

Previous work E! El El Q El El has shown that 
the state |NOPA) violates instances of the CHSH inequal- 
ity with the maximum violation being 2V2 El El El 
The violations in EH were arrived at by first applying 
displacement operations to modes 1 and 2 and then mea- 
suring whether each contained an odd or even number 
of photons. As all odd (even) numbers are represented 
by binary strings for which the least significant bit is 
one (zero), E3' s results tell us that \?Pcm), which is ab- 
stractly the same as |NOPA), violates the CHSH inequal- 
ity Si < 2 for the least significant bits in the binary 
representations of J\f\ and A2 when we apply appropri- 
ate displacement operations and then measure these bits 
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using the scheme in Subsection IIII Bl This fact is high- 
lighted in Fig. [21 which is a plot of results formally equiv- 
alent to those in 0] for the state IV'cm)- In particular, 
for the displacements a — (3 = and a' = —f3' = J, 
where J € 3?, it is a graph of S± versus J for squeezing 
parameter values of r = 0.5, r = 1 and r — 1.5. 
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FIG. 2: Plot of Si (J) versus J for r = 0.5, r = 1 and r = 
1.5 with displacements. The horizontal dotted line represents 
S%(J) = 2. Both Si (J) and J are dimensionless. 



B. Second least significant bits 

In addition to violating the inequality Si < 2, the state 
IV'Cm) also simultaneously violates the bitwise Bell in- 
equality 52 < 2. This can be seen by calculating the 
average E [02(21)62(22)] for this state, which is 



E [a 2 (zi)b 2 (z 2 )] 



l-2(Pr(02 =+l,&2 = -l|*l,S2) 

+Pr{d 2 = -1,6 2 = +l|zi,z 2 )) , 

(14) 



where Pr(a 2 — f,b 2 = g\zi,z 2 ) is the probability that 
a 2 = f and b 2 — g given the displacements Di(zi) and 
D 2 {z 2 ). As the second least significant bit of a bit string 
is '0' for the decimal numbers 0, 1, 4, 5, 8, 9 . . . and '1' oth- 
erwise, 



n 2 \zi,z 2 ) 



E [02(21)62(22)] = 1- 

2[ £ Pr(Mi=m,M 2 = 

\m=0,l,4,5,... 713=2,3,6,7... 

+ Yl J2 Pr(M x = ni,M 2 =n 2 \zi,Z2) I , 

ni=2,3,6,7... n 2 =0,l,4,5... / 

(15) 

where Pr{Mi = ni,M% = n 2 \z x ,z 2 ) is the probability of 
observing \Mi = ni)i\M 2 = n 2 ) 2 given the displacements 
Di(zi) and D 2 {z 2 ). This is known to be [3(j 

Pr(Ni = ni,J\f 2 = n 2 \zi,z 2 ) 



tanh p r p\ 



,1/2 



coshr 5! 
x exp(-(z 1 /xi 



ni - 

Mi 



V2 2 



tanh r 



>W)/2)| 



(16) 



A*i = 21 + 
tanh r and Lp 9 ^ is a generalised 



where p = mm(ni,n 2 ), q = max(ni,7i2 
z 2 tanhr, ji 2 = z 2 + 2* 
Laguerre polynomial. 

Calculating S 2 using Eqs l|15|l and l|lrj|) we obtain, upon 
setting a = 5 = and on a' = —f3' = J, where J 6 5f, 
CHSH violations for a range of J values. These are illus- 
trated in Fig.|3|as a function of J for squeezing parameter 
values of r = 0.5, r = 1 and r = 1.5. As is the case for the 
graphs in Subsection IIV CI and Section 52 was calcu- 
lated using Mathematics with all numerical errors being 
negligible |3jj ■ Significantly, for a range of J and r values 
we simultaneously violate the bitwise Bell inequalities for 
the least and second least significant bits in the binary 
representations of Mi and Ml, as can be determined by 
inspecting Fig. and results in [l4|. 



1.5 





(a) 




r=0.5 




r=T P- 


r=1.5 





0.05 0.1 0.15 0.2 0.25 

J 



2.08 




0.05 



0.1 j 0.15 0.2 0.25 



FIG. 3: 

(a) Plot of S2(J) versus J for r = 0.5 (solid line), 
r = 1 (dashed line) and r = 1.5 (-•-•-.) with dis- 
placements. Both S2 (J) and J are dimensionless. 
The horizontal dotted line represents 52 (J) = 2. 

(b) Close-up plot of S2(J) versus J for r = 0.5 with 
displacements. The horizontal dotted line represents 
S 2 (J) = 2. 
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C. Third least significant bits 

To show that |i/'cm) violates the bitwise Bell inequality 
S3 < 2, we now perform a similar calculation to the last 
subsection's except that, as the third least significant bit 
of the numbers 0, 1, 2, 3, 8, 9, 10, 11 ... is '0', 

E[a 3 ( Zl )b 3 (z 2 )] =l-2x 

EEE I>(M 

, i=even j—0 /—odd s— 



4i+j,A^ 2 =4Z + s|zi,z 2 ) 



EE E E Pr W = 4!+ ^ 2 = 41 + s i zi ' Z2 ) 



i— odd j—0 /—even s—0 



(17) 



Using this result and Eq. (jlHfl to calculate S3 as a func- 
tion of J, we obtain the results in Figs0](a) and (b) which 
show bitwise Bell inequality violations. Observe that for 
certain J and r values, we can simultaneously violate the 
bitwise Bell inequalities Si < 2, S 2 < 2 and S 3 < 2. It is 
also interesting to note that the violations in Fig. 0] are 
significantly less than those for the second least signifi- 
cant bits shown in Fig. It is possible that this is due to 
the fact that groups of four consecutive numbers (eg. 0, 
1, 2 and 3) share the same value for their third least sig- 
nificant bits. (For the second least significant bits, only 
two consecutive numbers share the same value.) Because 
of this, it may be more difficult for the displacement op- 
erations we implement to cause states to 'flip' the values 
of their third least significant bits. In turn, this would 
mean that it would be more difficult for these operations 
to generate the sort of interference between previously 
orthogonal states in IV'cm) necessary for obtaining Bell 
violations, thus leading to the smaller violations for the 
third least significant bits shown in Fig 01 



V. BITWISE BELL VIOLATION WITH LOCAL 
SQUEEZING OPERATIONS 

In this subsection we show that unitaries other than 
displacement operations can yield a bitwise Bell inequal- 
ity violation for |?/>cm)- In particular, we show that 
squeezing operations applied to the centre-of-mass vibra- 
tional states of both sets of ions in the x direction can 
produce such a violation involving the second least signif- 
icant bits of N\ and A/2. These squeezing operations are 
interesting to consider as they have been practically im- 
plemented in ion traps (see, for example, 35]). Observe, 
however, that squeezing operations applied to both sets 
of ions do not produce CHSH inequality violations for the 
least significant bits in the the binary representations of 
N\ and A/2 as squeezing operations are associated with 
two-phonon creation and annihilation. Thus, they do 
not cause odd and even phonon number states to change 
parity and so do not induce the type of interference re- 
quired for such violations. Throughout this section we 



2.002 
2.001 

r 

1.999L 
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2.1 

^2 
CO 

1.9 
1.8 
1.7 



0; 



0.05 0.1 j 0.15 0.2 0.25 



(b) 



r=1 



r=1 .5 



0.05 



0.1 j 0.15 



0.2 0.25 



FIG. 4: 

(a) Plot of S^J) versus J for r = 0.5 with dis- 
placements. Both Sz(J) and J are dimensionless. 
The horizontal dotted line represents Sz(J) = 2. 

(b) Plot of S3(J) versus J for r = 1.5 (solid line) and 
r = 1 (dashed line) with displacements. The horizontal 
dotted line represents Ss{J) = 2. 



assume that N > 2, so that the measurement scheme in- 
volves enough two-level electronic systems to to measure 
a 2 and 62 • 

Applying the above mentioned squeezing operations to 
I^Cm), we obtain 



\TPs)=S 1 (r + )S 2 (r-)S 12 (r)\0,0), 



(18) 



where S\ and S 2 are single-mode squeezing opera- 
tors for the centre-of-mass modes of the first and sec- 
ond sets of ions in |^cm) m the x direction with 

and r_. The opera- 



real squeezing parameters 
tor Si(r+) = cxp (r+ (bQ f 



(i) 



whilst 



S 2 (r_)=exp r_ [B. 



>(i) t 



B. 



(1; 



We now de- 



termine fV sqU eeze(A/i = ni,A/2 = n 2 |r + ,r_), the proba- 
bility of observing ri\ and n 2 centre-of-mass phonons in 
the x direction for the first and second sets of ions re- 
spectively given S\{r+) and S%(rJ) by re-ordering the 
operators in Si (r-i_)52(r , -)5i 2 (>-)|0)i0)2. The idea for 
this derived from a calculation in [3(| that found (A/i = 
ni,Af 2 = n 2 \ipD) by decomposing and normally ordering 
the operators in Di{a)D 2 ((3)S 12 (r)\ti) l {)) 2 . 
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Decomposing S\ 2 in a normally-ordered manner and 
utilizing the known single-mode squeezing operator de- 
composition i^H 

S,(R) = exp[-ln(co S hR)(B$lB% + 1/2)] x 

exp [tanh P cosh 2 R (b^J ^ /2] x 

expftanhi^B^) 2 ^], (19) 
where j = 1,2, yields 



exp (- lncosh r+ B<£ ] - lncosh r _ B$ f b£> 



tanhr_ cosh 2 r_ ( 1 ] /2 ] x 



exp ^tanhr + (B^Jj /2 + tanh 

exp(tanhrB^ ) t )]0, 0>, 
where if = cosh 2 r cosh r + cosh r_ 



K0 2 /2)x 



(20) 



To determine Pi 



squeeze 



(A/i = n\,N 2 = n 2 |r+,r_) from 
Eq. 120(1 we use three operator identities. These identi- 
ties, which hold for any operators A and B such that 
[A, A*] = [B,5t] = 1 and [A, i?t] = = 0, are 

proven in Appendix A and are 

exp(ciA 2 )exp(c 2 A t B t ) = exp(c 2 A t B t ) exp(cic^B t 2 ) 
x exp(2cic 2 AB t ) exp(ciA 2 ), 

(21) 



exp(ciA 2 )exp(c 2 A t2 ) = exp( 



e 2 



1 - 4cic 2 



exp [cosh _1 (l + - 



2cic 2 
- 4cic 2 



2c lC2 ){AU+l/2}] 



exp( 



ci 



1 - 4cic 2 



A 2 ) 



(22) 



and 



exp(ciA t B)exp(c 2j B t2 ) =exp(c 2j B t2 ) exp(2cic 2 A t B t ) 

x exp(c 2 c 2 A t2 ) exp(ciA t B), 

(23) 

where c\ and c 2 are complex numbers. Re-ordering terms 
on the right-hand side of Eq. (|20jl by applying these iden- 
tities to b[ x J and B 2 X ] yields 



\4>s) 



(24) 



exp ^— lncoshr+sj^ ^ b[ x J — lncoshr_£>, 
exp [(—1/2 tanh r + cosh 2 r + + 1/2 tanhr_ tanh 2 



exp 



-1/2 tanh r_ cosh 2 r_ + d 2 ) (-B^ f )' 



exp ((tanhr + 2d 1 d 2 )B { 1 1 ] f B<£> f )j |0,0>, 

where d\ — tanhr tanh r_, d 2 = d^/l — Ad^d^, d% — 
1/2 tanh r_, d^ — 1/2 tanh r + tanh 2 r and 

M = exp [cosh" 1 (1 + 2d3 f 4 - 2d 3 d 4 )l. 

1 — d^d^ 

Calculating Pr squee ze(A/i = ni,J\f 2 = n 2 |r + ,r_) using 
the right-hand side of Eq. 124fl . we arrive at 

-Pr sq ueeze(A/'i = ni,N 2 = n 2 |r + ,r_) = /(ni,n 2 ) x 



1 M 

— vni!n 2 !cosh ni r + cosh " 2 
if 



ni/2 n 2 /2 

r_e, e 



1 



min(ni,n2) 

,=o + /^4 + /... ( V^ / - j)/2)!((na - i)/2)l 



where 



/(ni,n 2 ) 



, rti + ri 2 = even 

1 , n\ + rt 2 = odd, 



(25) 
(26) 



tanh 2 



r + d 2 d 2 , 
tanhr + 



e + = —1/2 tanh r + cosh 2 r + + 1/2 tanh? 
e_ = — 1/2 tanh r_ cosh 2 r_ + d% and et v 
2d\d 2 . Observe that it was crucial to re-order the op- 
erators on the right-hand side of Eq. (|20fl in arriv- 
ing at Eq. I|25(l as if we did not then we would have 
had to deal with infinitely many terms contributing to 
fsqueeze(A/i = n\,N% — n 2 |r + , r_). The reason for this is 
that the right-hand side of Eq. I|20l) contains annihilation 
operators to the left of creation operators for the same 
mode. As a consequence, for example, upon consider- 
ing power series expansions of the exponentials in this 
equation we have contributions to |0,0) from terms in 
which we first create X centre-of-mass phonons, where 

X = 1,2,3..., by applying f-B^?^ , and then an- 



nihilate them by applying (Mx)^ to ( B i{x)) I ' )- 
Given that X can be any natural number, it follows 
that, to determine Pr sque eze(A /r i = 0,A/" 2 = 0|r + ,r_) us- 
ing Eq. H20f) . we seem to need to consider infinitely many 
contributing terms. In contrast, the only annihilation 
operators that appear to the left of creation operators 
in the right-hand side of Eq. I(24|l are present in terms 
containing number operators. These do not increase or 
decrease the number of centre-of-mass phonons when ap- 
plied to a state and so their presence does not cause 
infinitely many terms to contribute to P S qucozc(A/i = 
fiij-A/ja = n 2 |r + ,r_), thus making the calculation of 
P S queeze(M = Jiij-A/a = n 2 \r + ,r^) tractable. Using 
Pr squeeze (A/i = ni,N 2 = " 2 |r + ,r_) to calculate CHSH 
correlations in a similar manner to that used to deter- 



x 



mine Pr(N\ — rii,Af 2 — n 2 \zi,z 2 ) in Subsection IIV Bl 
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FIG. 5: 

(a) Plot of 52 (J) versus J for r — 0.5 with local squeezing 
operations. Both S2(J) and J are dimensionless. The hori- 
zontal dotted line represents the equation 52 (J) = 2. 

(b) Plot of £2 (J) versus J for r = 1 (solid line) and r = 1.25 
(dashed line) with local squeezing operations. The horizontal 
dotted line represents the equation 5a (J) = 2. 

we find that the CHSH inequality is violated, as illus- 
trated in Figs (a) and (b) which show £2 as a function 
of J for r — 0.5, r = 1 and r = 1.25. 



VI. DISCUSSION 

Though this paper's results are related to those in 
[3, they differ from 0's results in a number of ways. 
First, motivated by the comment in ^| that the formula- 
tion of a more practical measurement scheme to measure 
[ill's Bell inequalities was desirable, we arguably pro- 
posed such a scheme (at least for the least, second least 
and third least significant bits of Afi and A/2) which cen- 
tred on transferring the centre-of-mass vibrational state 
of group of trapped ions to the internal states of a group 
of electrons. Second, we violated inequalities similar to 
those in IIJ using vastly different schemes to those in 
|13| . In |13| the measurements made in violating the in- 
equalities were pseudo-spin measurements along varying 
axes in a two-dimensional plane. In contrast, we mea- 
sured a pseudo-spin based observable in a single direction 
and obtain violations by applying a range of displace- 



ments and squeezing operations to sets of N ions. 

Aside from extending work in the Bell violations 
in Section Hvl are noteworthy as they are similar to those 
attainable in so-called hyper- entangled states [39j . These 
are states in which more than four degrees of freedom 
are entangled, such as a two-photon state with polarisa- 
tion, energy and momentum entanglement (each particle 
has a polarisation, an energy and a momentum degree of 
freedom participating in the entanglement). As a conse- 
quence, hyper- entangled states can violate multiple Bell 
inequalities involving, collectively, five or more degrees of 
freedom. The Bell violations in Section Hvl are similar to 
those achievable in hyper-entangled states as the viola- 
tions in Section llVl involve violating three Bell inequal- 
ities involving six degrees of freedom, namely the three 
least significant bits of both A/"i and A/2 • One reason why 
this connection is interesting is that indirect evidence 
suggests [23 (pp. 2179-81) that some hyper- entangled 
states may be able to perform certain interesting quan- 
tum information processing. This, in turn, suggests that 
IV'cm) may also be able to perform such feats. 

Though we only demonstrated bitwise Bell violations 
for the three least significant bits of A/i and A/2 oth- 
ers, presumably, also exist for the fourth, fifth, sixth 
etc. least significant bits. However, the calculations 
required to demonstrate these violations were not per- 
formed as calculating S y becomes increasingly difficult as 
y increases due to the presence of more and more com- 
plicated spreads of centre-of-mass number states sharing 
the same value for the y th bit. An example of this in- 
creased complication can be seen by observing the fact 
that Eq. JT^) (y = 2) is simpler than Eq. JUJ (y = 3). 

An alternate approach we could have taken to investi- 
gating \ipcuYs nonclassical correlations would have been 
to see what quantum information processing tasks this 
state's correlations could be used to perform. However, 
one complication with this is that the 2N systems in the 
total physical system described by \tPcm) are not qubits 
but instead are infinite-dimensional harmonic oscillators. 
Because of this, we cannot directly consider if this state 
is useful in helping to implement well-known quantum 
protocols for qubits. In spite of this difficulty, however, 
a recent result showing that certain Bell violations im- 
ply the existence of quantum communication complexity 
protocols superior to any classical ones [2(| may be useful 
in manifesting |V>Cm)'s nonclassical correlations. In par- 
ticular, it may allow us to readily show that Section llVl s 
Bell violations imply that IV'Cm) could be employed to 
perform such quantum protocols. 

Yet another approach that could be taken to illustrate 
the nonclassicality of |?/>cm) is to use entanglement wit- 
nesses )40| . An entanglement witness W for the entan- 
gled state p is an operator such that Tr(pW) < and 
Tr(aW) > whenever a is a separable state. This ap- 
proach would involve identifying suitable operators W 
and then applying them to IV'Cm)- We acknowledge that 
it may be a useful approach to try, however, we have not 
explored it. 
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How feasible are the system and measurement scheme 
we have discussed? To reiterate, first, the system involv- 
ing the parametric oscillators feeding into two cavities 
within which lie ion traps containing N ions does not 
seem to be infeasible. This is because, as stated ear- 
lier, optical cavities and nondegenerate optical paramet- 
ric amplifiers have been widely realized in laboratories 
for some time. In addition, experiments in which a sin- 
gle harmonically trapped ion has been placed within an 
optical cavity have been conducted (27J- Another factor 
consistent with the potential feasibility of the system con- 
sidered is that the entangled state (i/'cm) can be created, 
to a good approximation, on a timescale far shorter than 
that of the vibrational decoherence for the ions. Finally, 
displacement and squeezing operations on trapped ions 
have been realised in [35j | via shining laser beams on the 
ions. 

To conclude, following on from 01 we have presented 
Bell inequalities that reveal certain nonclassical correla- 
tions in IV'cm)- In particular, these correlations are be- 
tween bits in the binary representations of M\ and A2 
that violate three bitwise inequalites. We have also pre- 
sented a bitwise Bell violation for |^cm) involving local 
squeezing operations. 



Acknowledgements 

Both authors wish to thank Professor S. Braunstein 
for stimulating discussions. DTP thanks Professor T. 
Bracken and Dr J. Links for assistance with Lie algebras 
and also acknowledges the assistance of H. M. Wolffram. 



Appendix A: Operator identities using Lie algebras 



In this appendix we prove the identities 

exp(c 1 A 2 )exp(c 2 A^B^) = exp(c 2 yl t S t ) exp(cic^B t 2 ) 
x exp(2cic 2 AS t ) exp(ciA 2 ), 

(27) 

02 A^ 2 )x 



1 - 4C1C2 



exp(ciA 2 ) exp(c2^4' 2 ) = exp( 
exp [cosh -1 (l + — ^ 2 Cl c 2 ){A^A + 1/2}] 

1 — 4ciC2 



exp( 



" A 2 ) 



1 — 4ciC2 



(28) 



and 



exp(ci^ t B)exp(c 2j B t2 ) = cxp(c 2 B t2 ) exp(2cic 2 A t B t ) 

x exp(c 2 c 2 ^4 t2 ) exp(ciA t S), 

(29) 

where the c\ and c 2 are c-number co-efficients and A 
and B are bosonic centre-of-mass annihilation opera- 
tors for two non-interacting systems for which [A, A^] = 
[B, iJt] = 1. We do this using the Baker-Campbell- 
Hausdorff (BCH) formula |4l| 1 a technique that has been 



called the differential- equation approach 38, 42] and the 
fact that SU(1,1) has a two-dimensional matrix represen- 
tation. 

The BCH formula is 41] (p. 118) 

exp(„4)exp(,B) = (30) 
exp(A + B + l/2[A,B] + l/l2{[A,[A,B}] + [[A,B},B}} 
+l/48{ [B, [A, [B, A}}} + [[[A, B],A],B]} + ...), 

where A and B are abitrary operators. We use it in 
proving Eqs ((27(1 and 1)29(1 by employing it to convert 
their left-hand sides into a single exponential each. Next, 
we convert these single exponentials to the normally- 
ordered products of exponentials on the right-hand sides 
of Eqs (|2T| and (|2"9"|) using the differential-equation ap- 
proach, which we now explain. This approach disentan- 
gles or decomposes a single exponential with a sum in its 
exponent, into a product of a number of exponentials in 
the following manner: First, we multiply the exponent of 
the single exponential by a parameter t. Next, we equate 
the single exponential with this extra factor of t in its 
exponent to a product of exponentials for which each ex- 
ponent is some unknown function of t multiplied by a 
generator of a certain Lie group. The group is the same 
for all exponents and is one for which all terms in the 
exponent of the original single exponential are some con- 
stant multiplied by a generator of the group. In addition, 
each generator appears precisely once in the product of 
exponentials. To give an example, consider the single ex- 
ponential exp(ci^l t2 + c 2 A 2 ). Noting that A 2 and A^ 2 
are generators of SU(1,1), we multiply c\A^ 2 + c 2 A 2 by t 
and equate exp ((ciA^ 2 + c 2 A 2 )t) to the following prod- 
uct of exponentials: 

exp {h{t)A^ 2 ) exp (f 2 (t)(A^ A + 1/2)) exp (f 2 (t)A 2 ) . 

(31) 

Observe that in expression 1(31(1 that each exponent is 
the product of an unknown function and a generator of 
SU(1, 1). Furthermore, each of SU(1, l)'s generators ap- 
pears exactly once. Returning to the general case, finally, 
we calculate the unknown functions of t, and so complete 
the process of disentangling the single exponential, by dif- 
ferentiating both sides of the equation equating the single 
exponential with the extra factor of t in its exponent to 
the product of exponentials with respect to t, multiply- 
ing both sides from the right by the inverse of the single 
exponential and, lastly, equating operator coefficients on 
both sides. The paper |42( contains a detailed example 
in which the differential-equation approach is used. 

Using the BCH formula on the left-hand side of 
Eq. J23 yields 



exp(ciA 2 ) exp(c 2 ^4 t 5 t ) 

2 

= exp(ci^ 2 + c 2 A^B^ + c lC2 AB^ + ^-B* 2). 



(32) 



Noting that {A^ B\ A 2 , AB\B^ 2 } form a Lie algebra, we 
use the differential-equation approach on the right-hand 
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side of Eq. I|32[l to know that there exists a normally- 
ordered decomposition such that 



exp[(ciA 2 + c 2 A^B^ + Cl c 2 AB^ + ^B^ 2 )t] 

o 

= exp(/ 1 (t)AtBt)exp(/ 2 (t)Bt2) x 

exp(f 3 (t)B^A)exp(f 4 (t)A 2 ) 
= U, 



(33) 



where fx, fy, / 3 and f± are functions we now determine. 
Differentiating both sides of Eq. (|33|) with respect to t 
and then multiplying from the right by U~ l yields 

Cl A 2 + c 2 A*Bi + c lC2 AB^ + ^B^ 2 

6 

= f^B* + f 2 e^ AtBt Bt2 e -A>itflt 
+ h e flA " Bf e f2B ' 2 B^a e'^" e~ hA ' B ' 
+ U e flAi ' Bf e hB " 2 e /aSU B^ 2 e~^ B ' A 



(34) 



where fi denotes dfo/dt. Using the identity |43 (p. 162) 
cxp(A)£exp(-A) = B + [A.B] + 1/2\[A, [A, B]] + 

(35) 

where A and B are arbitrary operators, on the right-hand 
side of Eq. produces 



ciA 2 + c 2 A^B^ + c lC2 AB^ + 2 

6 

= htfBt + f 2 B^ 2 + f 3 (B^A - /iSt 2) 
~h(A 2 -2f 1 B^A + f 2 Bl 2 ). 



(36) 



Upon equating operator coefficients we arrive at four cou- 
pled differential equations for the f^s. Solving these and 
setting t=l yields 

2 

exp(ciA 2 + c 2 A^B^ + Cl c 2 AB^ + ^flt 2 ) (37) 

6 

exp(c 2 A t J B t ) exp(ci c%B ] 2 ) exp(2cic 2j B t A) exp(civl ). 
Recalling that 

2 

exp(ci A 2 + ca^tfit + c x c 2 AB^ + ^B^ 2 ) = (38) 
exp(ci A 2 ) exp(c 2 A t B t ) 

we arrive at Eq. (27)l. The identity in Eq. l{2"5|) can be 

obtained via a very similar calculation to that which we 
have just performed. 



To prove Eq. l(2"5|) , we first note that A 2 and A* 2 are 
generators of the Lie grou p SU(1,1). Related to this 
group, it is known that [3Sj 

exp(ciA )exp(c 3 (AtA + l/2)) exp^A^ 2 ) 
= exp(/3i ({ Cl })A^ 2 ) exp( / 3 2 ({c l })(A t A + 1/2)) 

exp(/3 3 ({ Cl })^ t2 ), (39) 
where /3i, /3a and /3 3 are, as yet, unknown functions. 
Noting that these functions are only determined by the 
commutation-relation structure of SU(l,l)'s generators, 
we follow HH and replace A 1 * 2 , (A 1 * A + 1/2) and A 2 by 
two-dimensional matrices with identical commutation re- 
lations. This leads us to making the following transfor- 
mations: At 2 -> 2cr + , A 2 -> -2cr_ and At A +1/2 -> <r 2 , 
where 



and 



1 






1 



1 

-1 



(40) 



(41) 



(42) 



Upon doing this, and setting c 3 = 0, the left-hand side 
of Eq. l(3lljl becomes 



.■XJ)(-2r,,T_).-xj)(2r_>rT + ) = { ^ M >! 1 



whilst the right-hand side transforms to 



exp(/3i2<7+) exp(/3 2 <7 z ) exp(-/? 3 2cr_) 

P-4f3 x f3 3 M 2(3 X M 
-2(3 3 M M 



(44) 



where P = cosh /3 2 + sinh j3 2 and M = cosh (3 2 — sinh /?2 ■ 
Equating matrix elements on the right-hand sides of 
Eqs |33) and J33J leads to 

Pi = c 2 /(l - 4cic 2 ) 

f3 2 = cosh -1 (1 + 2cic 2 /(l - 4cic 2 ) - 2cic 2 ) 
and 

/3 3 = Cl /(l-4ciC2) (45) 
and hence to Eq. (|28[l. 
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